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Abstract : Present quantum theory, which is statistical in nature, does not 
predict joint probability distribution of position and momentum because they 
are noncommuting. We propose a deterministic quantum theory which pre- 
dicts a joint probability distribution such that the separate probability distri- 
butions for position and momentum agree with usual quantum theory. Unlike 
the Wigner distribution the suggested distribution is positive definite. The 
theory predicts a correlation between position and momentum in individual 
events. 
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1. Introduction . Present quantum theory does not make definite predic- 
tion of the value of an observable in an individual observation except in an 
eigenstate of the observable. Application of quantum rules to two separated 
systems which interacted in the past together with a local reality principle 
(Einstein locality) led Einstein, Podolsky and Rosen 1 to conclude that quan- 
tum theory is incomplete. Bell 2 showed that previous proofs of impossibility 
of a theory more complete than quantum mechanics 3 (inappropriately called 
hidden variable theory) made unreasonable assumptions; he went on how- 
ever to prove 4 that a hidden variable theory agreeing with the statistical 
predictions of quantum theory cannot obey Einstein locality. 

Bell's research was influenced by the construction by De Broglie and 
Bohm 5 (dBB) of a hidden variable theory which reproduced the position 
probability density of quantum mechanics but violated Einstein locality for 
many particle systems. For a single particle moving in one dimension with 
Hamiltonian 

H = -Ti 2 /{2m)d 2 /dx 2 + U(x), (1) 

and wave function ip(x,t), de Broglie-Bohm proposed the complete descrip- 
tion of the state to be {\(t), \ip)}, where \{t) is the instantaneous position 
of the particle, and its momentum is 

p dBB (X,t) = mdX/dt = [Re ^*(-ih / dx) / {\^\ 2 )\ x=x . (2) 

In an ensemble the position density p(X,t) agrees with |^(A,t)| 2 for all time. 
However, Takabayasi 6 pointed out that the joint probability distribution for 
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position and momentum given by the theory 

PdBB(A,p,t) = \^(\,t)\ 2 5 (p - p dBB (\,t)) (3) 

does not yield the correct quantum mechanical expectation value of p n for 
integral n ^ 1. De Broglie 5 stated that these values in his theory "cor- 
respond to the unobservable probability distribution existing prior to any 
measurement" and measurement will reveal different values distributed ac- 
cording to standard statistical quantum mechanical formula. Such a central 
role for measurement is unsatisfactory if one wishes to apply the dBB theory 
to closed quantum systems. 

Without using hidden variables, Griffiths 7 and Gell-Mann and Hartle 8 
introduced joint probability distributions for noncommuting observables at 
different times in the consistent history approach to quantum theory of closed 
systems. Wigner 9 had earlier introduced a joint distribution for x and p at 
the same time, 

/°° dy ( y \ ( y \ 

which yielded the correct quantum probability distributions separately for 
x and p on integration over p and x respectively. The Wigner distribution 
cannot however be considered a probability distribution because it is not 
positive definite, as seen from the fact that the integral 

J dx dp PwA x iP)PwA x iP) = lOM)! 2 /^ 2 ^) 
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vanishes for two orthogonal states ip, 0. 

We wish now to propose a deterministic quantum theory of a closed sys- 
tem with the following properties. (We consider in this paper only 1 particle 
in 1 space dimension). 

(i) At each time, the particle has a definite position x and a definite momen- 
tum p. 

(ii) The system point in phase space has a Hamiltonian flow with a c-number 
causal Hamiltonian H c (x, p, ip(x, t),t) so that in an ensemble of mental copies 
of the system the phase space density p(x,p,t) obeys Liouville's theorem 

dp(x,p,t)/dt = 0. (5) 

Here i/)(x,t) is the solution of the usual Schrodinger equation 

ik dip(x,t)/dt = HiP(x,t) (6) 

with H being the standard quantum mechanical Hamiltonian for the system 
and He being determined from the following criteria. 

(iii) Each pure "causal state", i.e., a set of phase space points moving ac- 
cording to a single causal Hamiltonian He has phase space density of the 
deterministic form 

p(x,p,t) = \i>(x,t)\ 2 5(p-p(x,t)), (7) 

in which p — p(x, t) = not only determines p as a function of x, but also de- 
termines X clS db function of p at each time (step functions being allowed when 
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necessary). Eqn. (7) guarantees on integration over p the correct quantum 
probability distribution in x for any real function p(x,t). The function is 
determined from the requirement that on integration over x, p{x,p, t) should 
also yield the correct quantum probability distribution in p. That such a 
determination is possible and unique apart from a discrete 2-fold ambiguity 
will be a crucial part of the present theory. It is obvious that our p(x, t) will 
have to be different from the p(x, t) of de Broglie-Bohm theory, 
(iv) Since the quantum probability distributions for x and p in the statistics 
of many measurements are exactly reproduced, so are the standard uncer- 
tainty relations. However, the correlation between position and momentum 
in individual events given by p(x, t) is an additional testable prediction of 
the present theory. 

In Sees. II, III we describe the construction of the momentum p(x, t) 
and the causal Hamiltonian He, in Sec. IV applications to simple quantum 
systems, and in Sec. V conceptual features of the new mechanics. 

2. Construction of Joint Probability Distribution of position and 

momentum. We seek a positive definite distribution of the form (7) 
where p is a monotonic function of x 

e dp(x,t)/dx > 0, e = ±l (8) 

The monotonicity property ensures that for a given t, the (^-function es- 
tablishes a one-to-one invertible correspondence between x and p whenever 
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dp/dx is finite and non-zero. The requirement of reproducing the correct 
quantum probability distribution of p is that 

j°°j{x 1 p ) t)dx = ^ (|,f) | 2 , (9) 

where ip(k : t) is the Fourier transform of ip(x,t). We substitute the ansatz 
(7) into (9) and integrate over momentum to obtain 

f dp' [ dx'\i;(x',t)\ 2 6(p'-p(x',t))= f &$(* t )\ 2 . (10) 

j — oo jp(x' ,t)<p J -co n \n j 

The region p(x', t) < p becomes x' < x if e = 1, and x' > x if e = — 1, where 
p(x,t) = p. Thus, we obtain, for e = ±1, 

^'^(ex',^! 2 = / rfA;'|^(A; , ,t)| 2 . (11) 

The left-hand side is a monotonic function of x which tends to 1 for ex — > oo 
for a normalized wave function; the right-hand side is a monotonic function 
of p tending to 1 for p — > oo (Parseval's theorem). Hence, for each t, Eq. 
(11) determines two monotonic functions p of x, one for each sign of e. (Note 
that the curve p(x, t) may have segments parallel to x-axis or p-axis corre- 
sponding to i/)(x,t) or ip(p/h,t) vanishing in some segment). The two curves 
p = p±(x,t) so determined yield via Eq. (7) phase space densities p±, with 
different causal Hamiltonians {Hc)± determined below. 

3. Determination of the Causal Hamiltonian. We view p(x, p, t) as de- 
scribing an ensemble of system trajectories in the phase space. We saw in the 
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last section that such a description is possible at each time. We would now 
like to find a causal Hamiltonian such that the time evolution in phase space 
implied thereby is consistent with the time dependent Schrodinger equation. 

In order that the total number of trajectories is conserved in time we 
must have the continuity equation 

dp/dt + d(px)/dx + d(pp)/dp = (12) 

If the dynamics of the trajectories is of Hamiltonian nature i.e. 

x = dHc/dp, p = -dHc/dx (13) 
then we have Liouville's theorem that the phase space density is conserved, 

dp/dt + xdp/dx + pdp/dp = (14) 

i.e. 

dp/dt + (dHc/dp) dp/dx - (dH c /dx) dp/ dp = 0. (15) 

The c-number Hamiltonian H c describing the causal time evolution of the 
trajectories in the phase space will be allowed to be different from the usual 
g-number Hamiltonian H describing the time evolution of the Schrodinger 
wave function tp according to Eq. (6). 

On substituting into Eq. (15) the ansatz (7) discussed in the last section, 
we obtain 

£S(p-p) + ^(ri6(p-p))=0 (16) 
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where 

f = d\ip\ 2 /dt + (dHc/dp) d\ij\ 2 /dx - drj/dp 
Tj = — l^l 2 {dp/dt + (dp/dx) dH c /dp + dH c /dx} . 
We thus need for consistency 

£ = and 77 = if p = p. (17) 

We now specialise to the usual case when H is given by Eq. (1). We find 
that this situation is taken care of with the choice of Hc(x,p, t) 

H c = ^-(p-A(x,t)) 2 + V(x,t). (18) 

The causal Hamiltonian is of the Newtonian form apart from the introduction 
of a vector potential A(x,t) and allowing the potential V(x,t) to differ from 
U(x). Eqs. (17) lead to the following equations to determine V and A (after 
using Schrodinger eqn. to substitute for d\tp\ 2 /dt), 

-dV(x, t)/dx = dp(x, t)/dt + (2my 1 d (p(x, t) - A(x, t) f /dx, (19) 

d[\^\ 2 (p- A-mv)] /dx = 0, (20) 

where v is given by 

v(x,t) = h/(2im) d£n(il)/ij)*)/dx (21) 

which is just the de Broglie-Bohm velocity. Eq. (20) implies that the quantity 
in square brackets must be a function of t alone. We choose this function 
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of t to be zero in order to avoid a singularity of the vector potential at the 
nodes of the wave function. We thus obtain 



A(x, t) = p(x, t) — mv(x, t) 



(22) 



With the calculation of the causal Hamiltonian thus completed via Eqs. (18), 
(19) and (22) a consistent Liouville description emerges. 

4. Illustrative Examples, (i) Quantum Free Particle. Let the quantum 
free particle be described by the Gaussian momentum space wave function 

${p/%, t) = (2tt)- 1/4 exp [-(p - [3) 2 /(2ah 2 ) - ip 2 t/(2mh)] (23) 

so that the coordinate space wave function is 



4>(x, t) = (ira) (ma/ (m + iaht)) 1 ^ 2 exp /, 



(24) 



/ = -(«/2) 



[x - fit /my - i x 2 + 



aht 2 2px p 2 t \ 



m 



ah mah 



1 + 



a 2 h 2 t T 



Our procedure yields 

A = p - /3 = ±h } 

and 



m 2 a 2 



+ a 2 k 2 t 2 \ X ml' 



dV/dx = ±{m 2 + a 2 h 2 t 2 )- 3/2 [xt(ah) 2 + (im}(ham) 



(25) 



(26) 



The determination of the causal Hamiltonian is now complete apart from 
an irrelevant additive function of t. Eq. (25) apart from predicting the 
momentum p — (3 at the centre of the wave packet x = (3t/m, also predicts 
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values of p at other values of x which the particle may have in individual 
events. The quantum potentials A and V are seen to be proportional to Ti in 
this example. 

(ii) Quantum Oscillator. For the minimum uncertainty coherent state of 
the harmonic oscillator of mss to, angular frequency u and amplitude of 
oscillation a we find 



pffuj 



rmo 

t exp 
Tin 



Imu 2 
— (x — a cos(u;tJJ 
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S(p-p(x,t)), (27) 



where 



and 



p(x,t) = —mujasm(ut) ± muj{x — a cos(cjt)), (28) 



A(x,t) = ±mu(x — acos(cjt)), (29) 
—dV(x,t)/dx = — muo 2 a cos(ujt) ± muj 2 asm(ujt) (30) 
The causal Hamiltonian yields the equation of motion 

md 2 x/dt 2 = — muj 2 a cos out (31) 

which results in exact harmonic motion even for x away from the centre of 
the packet. We do not of course expect this for solutions of the Schrodinger 
eqn. different from the coherent state here considered. 

5. New conceptual features, (a) We have derived corresponding to every 
quantum wave function ip, two joint probability distributions for position and 
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momentum of the form (7) which are (i) positive definite, (ii) have Hamilto- 
nian evolution with causal Hamiltonians (H c )± and obey (hi) 



(f(x)+g(p))p±(x, P ,t)dxd P = UAf{x)+g\-ih—\ U , (32) 



for arbitrary functions f(x) and g(p). Eq. (32) is the major advantage of the 
present theory over the dBB theory. It raises the exciting possibility that 
the momentum values p±(x,t) for individual events here derived could agree 
with experimental values, and the single particle theory described could be 
the seed of a general quantum theory of closed systems. We postpone the 
discussion of measurements until we present a generalization of the theory 
to many particles. 

(b) Since both p + and p_ obey Eq. (32) so will p = C p + + (1 — C)p_ with 
< C < 1. But since p + and p_ correspond to different causal Hamiltonians 
(Hq)±, p will not correspond to a 'pure causal state'. We are led to the 
concept of a pure causal state as being more fine grained than a pure wave 
function ip. All p = Cp+ + (1 — C)p_ correspond to ip (p <-> ip) for a 
continuum of values of C, but only C — 0,1 correspond to pure causal 
states. To quantum density matrix states T:C a \tp a ) (tp a \ correspond phase 
space densities SC Q p Q if p a <-> ip a . 

(c) It is clear that the causal Hamiltonian evolution of phase space densities 
could be described purely in the phase space language without using the 
intermediate step of the wave function. We find it convenient to use ip(x,t) 
at the present stage. 
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(d) One can ask if Eq. (32) can be generalized to more general quantum ob- 
servables. Here we face the old problem that there exist nonclassical observ- 
ables e.g. x (— ih-^j x, ih-^j xx + h.c.^j /2 which have different expecta- 
tion values but the same 'naive' classical analogue x 2 p. A trivial way followed 
already for the dBB distribution is: given a nonclassical observable A the 
phase space analogue can be f(x,p,ip) such that f(x,p,ip) = ip*Aip/\ip\ 2 . 
Perhaps only those operators A, which like (f(x) + q{p)) have a phase space 
representation f(x,p) independent of ip should be considered as 'beables'. 10 

(e) The predictions of the momentum values p(x, t) are independent of any 
special ansatz for He- 

(f ) Due to the existence of trajectories, the problem of inconsistent histories 
will not arise in this theory. 

John Bell's writings, especially his last article 'against measurement' 10 
has profoundly influenced this work. We are grateful to colleagues at the 
Tata Institute especially to Deepak Dhar for many stimulating questions. 
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